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ABSTRACT 

Frequency glitches have been observed so far only in radio pulsars and anomalous X-ray pulsars. 
Here we present evidence of a glitch in a neutron star accreting from a Be companion. The transient 
KS 1947+300 reappeared in 2000 October as a moderately strong X-ray source that exhibited 18.7 s 
pulsations, leading to an identification with the BATSE source GRO J1948+32, last detected in 1994. 
We have analyzed Rossi X-ray Timing Explorer observations taken during the 2000-01 outburst, as 
well as additional observations taken during a smaller outburst in July 2002. Orbital Doppler shifts 
are apparent in the temporal variation of the pulse frequency. A joint fit of the RXTE data together 
with BATSE measurements from an outburst in 1994 yields the orbital period P or t, — 40.415±0.010 d, 
the projected orbital radius a x sin i = 137 ± 3 lt-s, and the eccentricity e = 0.033 ± 0.013. This degree 
of eccentricity is unexpectedly low for such a wide orbit. Pulse timing results also show that the 
intrinsic pulse frequency increased from 53.30 to 53.47 mHz at a rate approximately proportional to 
the X-ray flux. This is about the degree of spin up expected from the accretion torques that must 
be present when the X-ray luminosity reaches ~ 10 38 ergs s — 1 . On one occasion during the 2000-01 
outburst, the pulse frequency increased by ~ 1.8 x 10 -6 Hz in < 10 hr over and above the mean trend 
seen around that time, without any indication of a correspondingly large increase in X-ray flux. The 
fractional change in frequency of 3.7 x 10~ 5 during this event is significantly larger than the values 
observed in the glitches in radio pulsars and anomalous X-ray pulsars. We discuss other similarities 
and differences between these events. 

Subject headings: accretion — pulsars: general — pulsars: individual (KS 1947+300) — X-rays: stars 



1. INTRODUCTION 

Abrupt changes ( "glitches" ) in the spin frequency of ro- 
tating neutr on stars have been ob served to date in radio 
pulsars fe.g. lKrawczyk et alJl200 3fl and anomalous X-ray 
pulsa rs (AXPs: e.g. iKaspi et al.l l2003t IKaspi fe Gavriill 
2003). These events appear as sudden increases in the 
observed pulse frequency, superimposed on more-or-less 
steady spin-down arising from magncto-dipole braking, 
and are attributed to a sudden transfer of angular mo- 
ment umlDetweejitl^stnierfh^ and the 
crust l|Rudermanll976t [Anderson & Itoh 1975). Between 
glitches the interior supcrfluid and the crust are relatively 
loosely coupled, and the braking torque acts primarily on 
the crust, so that just prior to a glitch the crust is spin- 
ning slower than the interior superfluid. 

To date there have been no reports of glitches in accret- 
ing pulsars. Here we present the results of a pulse tim- 
ing analysis of data obtained from Rossi X-ray Timing 
Explorer (RXTE) observations of the accreting pulsar 
KS 1947+300, from which we derive orbital parameters 
as well as the intrinsic spin frequency history of the neu- 
tron star. The latter includes a short time interval during 
which the pulse frequency increased at an unusually high 
rate, that may be the first evidence for a frequency glitch 
similar to those observed in radio pulsars and AXPs. 

The cosmic X-ray source KS 1947+300 (I 11 = 66?09, 
h n = +2 °io) was first detected on 1989 June 8 with 
the TTM coded mask i maging spectrometer mounted 
on M ir's Kvant module (jSkinner| Il989t iBorozdin et al.l 
1990). The optical counterpart is a V = 14.2 BOVe star 
with moderate reddening that indicates the distance to 
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the system is « 10 kpc l|Negueruela et all2003j) . In 1994 
the Burst and Transient Source Experiment (BATSE) 
aboard the Compton Gamma Ray Observatory (CGRO) 
detected 18.7 s pulsations from an X-ray source within 
a few degrees of KS 1947 +300. From the BATSE data, 
iChakrabartv et all 1(19951) found that the apparent pulse 
period of GRO J1948+32, as the pulsating source was 
designated, varied systematically over the 33 d outburst, 
indicating that the source was in a binary system with 
an orbital period between 35 and 44 days. In 2000 
November, weak emission from KS 1947+300 was de- 
tected with the All-Sky Monitor (ASM) aboard RXTE. 
A search for periodicities in the ~ 5 yr ensemble of 
ASM measurements of KS 1947+300 accumulated up un- 
til that time yielded evidence o f a 41.7 + 0.1 day peri- 
odicity l)Levine fc Corbeill200 0L Soon afterwards, pul- 
sations with a period of 18.7 s were found in follow- 
up observations with the Proportional Counter Array 
(PCA) on RXTE, which showed that KS 1947+300 and 
GRO J1948+32 are one and the same ((Swank fc Morganl 
12000(1 . A series of RXTE target-of-opportunity (TOO) 
observations were then conducted which, together with 
the measurements by the ASM, gave excellent coverage 
over a long and strong outburst. 

2. OBSERVATIONS AND ANALYSIS 

We made observations of KS 1947+300 with all 
three instruments aboard RXTE: the ASM, PCA and 
the High-Energy X-ray Timing Experiment (HEXTE). 
The ASM views m uch of the sky every few hours 
ijLevine et al.l li996). The instrument consists of three 
Scanning Shadow Cameras (SSCs) mounted on a rotating 
platform, and carries out sequences of 90 s observations 
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Fig. 1.— Intensity of KS 1947+300 throughout the 2000-01 out- 
burst, (upper panel) The gray markers show the ASM count rate 
(1.5-12 keV) la confidence regions (with scale on the left in ASM 
counts s — 1 ), from 1 day averages of 90 s dwells, smoothed over a 
5-d window (excluding those with estimated Ict errors greater than 
0.7 count s -1 )- The open squares show the integrated 2—60 keV 
source flux (with scale on the left in units of 10 -9 ergs cm -2 s _1 ) 
as measured from best-fit models of PCA and HEXTE spectra. 
MJD 52000 corresponds to 2001 April 1. (lower panel) The flux 
measurements are plotted as squares on an expanded time axis, 
along with the individual Gaussian components (dotted lines) used 
to m odel the flux evolution, and the total model (dashed line; see 
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called dwells, in between which the platform holding the 
SSCs is rotated by 6°. The data from each SSC are 
analyzed to obtain estimates of the intensities of all the 
sources listed in the ASM analysis catalog within the field 
of view. The intensities are expressed as the count rate 
expected from a source at the center of the field of view 
of SSC 1 in 1996 March 1 

The ASM 1.5-12 keV intensity measurements of 
KS 1947+300 from 2000 July through 2002 April are 
shown in Fig. ^ The source became active in 2000 Octo- 
ber and rose to an intensity of approximately 1 count s -1 
around 2000 November 7 (MJD 51855). Following this 
initial peak the intensity declined somewhat before be- 
ginning a second, more substantial increase which peaked 
around MJD 51950 at ~ 9 counts -1 (approximately 

1 The Crab Nebula appears in the ASM as a 75 count s — 1 source 
(1.5-12 keV) and has a corresponding energy flux in that energy 
band of ~ 2.7 X 10 -8 ergs cm~ 2 s -1 . 



120 mCrab). Following the outburst maximum, the 
source declined to < 1 count s -1 , but subsequently exhib- 
ited a number of small outbursts peaking at 1-3 count s -1 
and lasting for 10-50 d. The source was not detected 
by the ASM prior to MJD 51800; the intensity was 
< 0.2 counts -1 . Relatively few measurements were 
taken around MJD 51960, and again around MJD 52290, 
when the Sun was relatively near its closest angular dis- 
tance of ~ 50° from the source. In addition, the mea- 
surements taken close to those times were more strongly 
affected by systematic problems than measurements at 
other times. 

In 2000 November, shortly after the outburst was no- 
ticed in the ASM data, a series of observations with the 
PCA and HEXTE aboard RXTE was initiated. The 
PCA consists of 5 Proportional Counter Units (PCUs) 
sensitive to X-ray photons in the range 2.5-90 keV, each 
with a collecting area of ~ 1400 cm 2 and a 1° field 
of view l)Jahoda et al.l[l996l) . Photon arrival times are 
measured to ~ 1 fis. Spectra are accumulated in up 
to 256 energy channels. The number of PCUs oper- 
ating during the observations varied between 2 and 5. 
The HEXTE comprises two clusters, each of which con- 
tains 4 scintillation detectors sensitive to photons in the 
range 15-250 keV, c ollimated to view a common 1° field 
l)Cruber et al.lll996T) . The detectors in the two clusters 
provide a total collecting area of 1600 cm 2 . 

The first PCA/HEXTE observation of KS 1947+300 
was performed on 2000 November 21 (MJD 51869). Sub- 
sequently, 1-6 ks long PCA/HEXTE observations were 
made every 1-2 days until the end of the main out- 
burst around 2001 June 18 (MJD 52078). Near the peak 
of the main outburst two long observations were made 
(for proposal ID 50068, PI: W. Heindl) that yielded 51 
and 98 ks of exposure, respectively. In order to obtain 
a longer baseline for determination of the orbital pe- 
riod, an additional series of observations began on 2002 
March 29 (MJD 52362) and continued until 2002 May 17 
(MJD 52411). In all, 132 observations were made with a 
total observing time of 371 ks. 

Spectral analysis of the PCA and HEXTE data was 
undertaken using lheasoft release 5.1 (2001 June 26). 
The data were first screened to ensure that the center of 
the fields of view were within 0.02° of the direction of 
KS 1947+300 and that the limb of the Earth was more 
than 10° from the source direction. Mean energy, i.e., 
pulse height, spectra were extracted for each observa- 
tion from standard observing modes ( "Standard-2" , with 
129 channels between 2-60 keV for PCA, and "Archive" , 
with 64 channels between 15-250 keV for HEXTE). PCA 
spectra were accumulated separately for each PCU. In- 
strument response matrices were generated for each PCU 
and each observation using pcarsp. We estimated back- 
ground count rates in the PCA using either faint or bright 
source models (depending upon whether the count rate 
exceeded « 40 counts s -1 PCU -1 ) developed for gain 
epoch 5 (from 2000 May 13 onwards) with pcabackest 
version 2.1e. We measured fluxes by integrating model 
fits to PCA and HEXTE spectra (see ©. 

Pulse timing analyses were performed using software 
written by the RXTE team at MIT. For timing anal- 
ysis, the 3-12 keV energy range was found to provide 
nearly optimal signal-to-noise ratios during the faintest 
observations in which we could detect pulsations. Thus 
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we selected, from "GoodXenon" mode data, only those 
events within the energy range 3-12 keV detected in the 
top layer of any PCU, adjusted the times of the events ac- 
cording to the projected distance of the spacecraft from 
the Solar System barycenter, and binned the selected 
events into 0.25 s time bins. 

3. SPECTRUM 

We fitted various continuum models to PCA spectra 
in the range 2-25 keV and HEXTE spectra in the range 
15-80 keV using xspec version 11.0. Above a flux of 
w 8 x 10~ 10 ergs cm -2 s _1 (2-60 keV) the best fit was ob- 
taine d with a combination of Comptonization ijTitarchukl 
119941 'compTT' in xspec) and blackbody components, 
along with a Gaussian representing fluorescent Fe Ka 
emission centered around 6.5 keV. Other models typi- 
cally used for pulsar spectra, including cutoff or broken 
power laws, or a power law with a blackbody component, 
generally gave worse fits. At high count rates, system- 
atic errors in the PCA response are known to contribute 
to deviations from the model; for example, the reduced- 
X 2 (xi) f° r the Comptonization model fit to the highest 
flux spectrum was 3.41 (i.e. x 2 — 511.7 for 150 degrees 
of freedom). Adopting a systematic error of 1.7% (2% is 
often used for the PCA) was sufficient to reduce the fit 
statistic to 1.01. For comparison, a fit with this assumed 
systematic error and the same Gaussian component com- 
bined with a power law with a high-energy exponential 
cutoff gave xl — 1-52, while a combination of blackbody 
and power law components gave 5.29. 

The PCA spectrum was significantly softer near the 
peak of the outburst compared to spectra obtained dur- 
ing the rise and decay (Fig. This was primarily due to 
a decreased contribution from the Tbb = 3-4 keV black- 
body component at higher fluxes, rather than a system- 
atic variation in Tbb (for example). While Tbb is unusu- 
ally high for a neutron star, examination of the unfolded 
spectra suggests that the requirement for the additional 
component may be due to intrinsic curvature of the ap- 
proximately power law part of the Comptonized spec- 
tra. A similarly good fit was obtained with a very broad 
(<7 « 10 keV) Gaussian component in place of the black- 
body, and we note that the integrated flux is not sensitive 
to the nature of the additional component. Regardless of 
whether a blackbody or Gaussian component is used, the 
additional component is strongly required for the best fit. 

The fitted column densities of neutral (unionized) 
matter measured over the outburst were consistent 
with the line-of-sight value est imated from Hi /dust sur- 
veys, i.e., ~ 10 22 ato ms cm -2 (jDickev fc Lockmai]ll990t 
iSchleeel et al.lll998|) . Thus, to obtain flux estimates we 
fixed Tin at this value and integrated the best-fit spectral 
models comprising Comptonization and blackbody com- 
ponents over the 2-60 keV band. The results are shown 
in Fig. ^ At the peak of the outburst in 2001 February 
the 2-60 keV flux was ~ 9 x 10 -9 ergscm~ 2 s _1 . The 
X-ray luminosity was therefore ~ 1.1 x 10 38 ergss -1 for 
a source distance of 10 kpc. 

4. PULSE TIMING ANALYSIS 

4.1. Frequency measurements 

To measure the long term frequency evolution of the 
pulsar we grouped together those RXTE observations 



3 




10 100 
Energy (keV) 

Fig. 2. — Representative photon spectra of KS 1947+300 from 
the peak (2001 February 10, or MJD 51950) and late in the decay 
(2001 May 3, MJD 52032) of the 2000-01 outburst. The effect of 
the instrumental responses have been removed ("unfolded"). Data 
from PCU #2 are shown, along with data from HEXTE clusters 1 
and 2. For the 2001 May 3 observation, we rebinned the HEXTE 
spectra by a factor of 4. Error bars indicate the lcr uncertainties. 
The top and bottom panels show the residuals to the model fit 
(units of normalized counts s _1 keV - 1 ) for the 2001 February 10 
and 2001 May 3 spectra, respectively. 

that began within 8 hr of each other and would form a 
combined data set spanning less than 24 hr. The fre- 
quency was then estimated for each set of observations 
by fitting the 0.25-s time series with the function 

F(t) = a +ai sin(wt)+a2 cos(wt)+a 3 sin(2(jjt)+a4 cos(2cjt) 

over a grid of trial frequencies where the coefficients ai 
were chosen to obtain the best fit for each frequency lu. 
The accuracy of each frequency measurement was esti- 
mated by the range of frequencies where x 2 < Xmin + 10 
(equivalent to ss 3<r), where x 2 m in was the value of the 
goodness of fit statistic for the best fit. 

We also used pulse frequencies determined from obser- 
vations performed in 1994 with the BATSE instrument 
on the Compton Gamma Ray Observatory. The frequen- 
cies were obtained i n the same way as those reported by 
iChakrabartv et aD l)1995|) except for being determined 
from 1-day rather than 2-day data sets. Data for the 
photon energy range 20-60 keV were used in this analy- 
sis. 

The measured frequencies are shown in Fig. |21 In 
the w 6 year interval between the 1994 and 2000 out- 
bursts the neutron star rotational frequency decreased 
from 53.47 to 53.30 mHz, corresponding to a mean spin- 
down rate / « -8.2 x 10~ 13 Hzs^ 1 . During the 2000- 
2001 outburst, there was a substantial overall increase 
in pulse frequency likely caused by accretion torques on 
the neutron star. An approximately sinusoidal periodic 
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Fig. 3. — Barycenter-corrected pulse frequency for KS 1947+300 
throug hout the 1994 Ap ril outburst as measured by BATSE 
IChakrabartv et al.|[l995l) . and the 2000-01 and 2002 outbursts 
as seen by RXTE. The upper panel shows the frequ ency measure- 
ments and the best-fit model (solid line; see H4.ll , The sum of 
model components representing the intrinsic neutron star spin fre- 
quency is shown as a dashed line. The lower panel shows the fit 
residuals. Error bars represent lc uncertainties. The grey his- 
togram shows the PCA/HEXTE flux measurements. 

variation is superposed on this increase; we naturally at- 
tribute it to Doppler shifts from orbital motion in a bi- 
nary system. 

4.2. Frequency Model 

We fit the frequency measurements with a simple 
model comprising the effects of neutron star spin includ- 
ing the spin-up due to accretion torques and the apparent 
changes due to orbital Doppler shifts: 

27r/oi a xSini 

/(*) = fspin(t) £ (cos / + g sin 21 + hcos2l) 

+ ( orb 

where f sp in(t) is the time-dependent stellar spin fre- 
quency, foj is an initial spin frequency for each of the 
three observation epochs (BATSE, RXTE 2000-01, and 
RXTE 2002) denoted symbolically by j, a x sini is the 
projected orbital semimajor axis in units of light travel 
time, and P rb is the orbital period. The coefficients 
g (= esinw) and h (— ecosw) are to be interpreted as 
functions of the eccentricity e and the longitude of pe- 
riastron u>. Finally, I — 2Tr(t — T^/ 2 )/ P m h + tt/2 is the 
mean longitude, with T„/ 2 the epoch at which the mean 
longitude I = tt/2. For a circular orbit, T w / 2 lS the epoch 
of superior conjunction of the neutron star. The right- 
most term in parentheses of Eqn. 2 represents the orbital 
Doppler shifts to first order in e; given the magnitudes 
of the uncertainties of our measurements, this should be 
an adequate approximation as long as e < 0.2. 

The neutron star spin frequency may be affected by 
several physical mechanisms. Among these, accretion 
torques are likely to be most important. As a first ap- 
proximation, we take the accretion torque to be pro- 



portional to a power of the observed X-ray flux, i.e. 
r acc oc Fj£. This prescrip tion is a limited generali za- 
tion of models like those of iGhosh fc Lambl l)1979altf) in 
which the torque is expected to be proportional to the 
accretion rate and, in turn, the flux raised to the power 
6/7. 

To obtain a simple model of the X-ray flux Fx, we 
fitted a composite light curve consisting of daily aver- 
ages of PCA/HEXTE 2-60 keV fluxes and ASM 1.5-12 
keV intensities (renormalized to 2-60 keV fluxes) with 
the sum of a set of Gaussian components. The 2000- 
2001 outburst was modelled by the sum of 5 compo- 
nents, while the time interval of the 2002 observations 
was modelled by a single component. The parameter 
values defining the model are listed in Tabled The in- 
dividual components that describe the 2000-2001 out- 
burst, as well as the total model, are plotted with the 
measured flux values in Fig. ^ lower panel. Over the 
2000-2001 outburst, the root-mean-squarc difference be- 
tween the measured fluxes and the model comprising 
five Gaussian components (nos. 2-6 of Table 1) was 
0.32 x 10 -9 ergs cm -2 s -1 . We also adopted an additional 
Gaussian component to represent the flux evolution dur- 
ing the 1994 BATSE observations. For that outburst, 
which lasted around 30 d, the source reached a maximum 
pulsed 20-60 keV flux of around 4.5 x 10~ 10 ergs cm -2 s _1 
on MJD 49465 l|Chakrabartv et alJll995l) . Broadband 
PCA/HEXTE spectra from the 2000-2001 outburst in- 
dicate that this likely corresponds to a 2-60 keV flux of 
around 1.9 x 10~ 9 crgscm _2 s _1 for a pulse fraction of 
50% (see i)4.3fl . Thus, we included a Gaussian compo- 
nent centered on MJD 49465 with a la width of 10 d 
and peak flux of 1.9 x 10~ 9 ergs cm -2 s _1 . 

We first tried to fit the measured frequencies with a 
simple model based on Eqn. 2 in which the time inte- 
gral of r acc gives the change in intrinsic spin frequency. 
The best fit was poor (\ 2 — 2715 for 101 degrees of free- 
dom). There are a number of possible reasons for this, 
including (1) the spacing of the observations does not al- 
low us to construct a flux model that reliably follows flux 
variations on short time scales (< 10 days), (2) the accre- 
tion rate may not be precisely proportional to the X-ray 
flux, and (3) the accretion torque in this simple model 
vanishes when Fx = 0. The latter property contradicts 
the observation of a net spin down during the w 6 years 
between the BATSE and RXTE observations while the 
source was in quiescence, i.e. while Fx ~ 0. In order to 
allow the net model torque to be negative when Fx = 0, 
we have augmented the spin model with terms linear in 
time separately for each of the three epochs. The revised 
model of the intrinsic spin frequency is thus: 

fs P in(t)=fo j + b j (t-t 0j )+l3[ F£(u)du+Af el (t) (3) 

Jt oj 

where foj and toj are the frequency and time, respec- 
tively, of the first frequency measurement for epoch j and 
Fx(t) is the 2-60 keV flux in units of 10 -9 ergscm _2 s _1 . 
The adjustable parameters foj and bj are set indepen- 
dently for each epoch while the power law index a 
and torque coefficient (3, along with the orbital param- 
eters, are adjustable parameters that apply globally to 
all three epochs. Further analysis (see i!4.3[) indicated 
that there was a rapid change in the intrinsic pulse 
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Fig. 4. — Representative pulse profiles from KS 1947+300 ranked 
by mean count rate. From top to bottom, the profiles are for 
the observations of 2001 Feb 10, Apr 5, May 6, May 12, May 17, 
May 18, May 21, May 27, 2000 Nov 28, and 2001 Jun 18 (MJD 
51950, 52004, 52035, 52041, 52046, 52047, 52050, 52056, 51876 and 
52078), respectively. The profiles have been rescaled by the mean 
rate and shifted in phase so that the minimum falls at phase 0.0. 
Each profile is plotted twice and shifted vertically for clarity. 

frequency around tgi ~ MJD 51956.41. We therefore in- 
cluded the term Af g \(t) which was assigned a value of 
zero for all measurements made prior to t g \ and a con- 
stant value of 1.83 /iHz for all measurements made af- 
terwards. The magnitude of this frequency change was 
determined through a pulse arrival time analysis which 
is described below, and so was not a free parameter in 
the frequency model. 

The model can be formulated so that it depends lin- 
early on all the adjustable parameters except for a and 
Port,. Thus, we searched a grid of points in a-P rb 
space to determine the best-fitting parameter combina- 
tion. When obtaining the best fit for each a-P rb pair 
of values, all other parameters were allowed to vary in a 
linear \ 2 fitting process. The values of the parameters 
for the best-fit model are listed in Tabled The model 
fits the data much better than our original simple model 
described above, although there are still significant resid- 
uals (Fig. This is reflected in the reduced x 2 value of 
3.35. The fit results unambiguously yield the number of 
orbital cycles from the epoch of the BATSE observations 
to that of the RXTE observations, and thus an accurate 
orbital period of 40.415 d. We also note that the values of 
the coefficients bj, i.e., the coefficients of the terms rep- 
resenting the values of / at zero flux for each epoch, are 
negative as expected from the spin down inferred during 



quiescence. 

We first estimate the uncertainty in each fitted param- 
eter of the frequency model from the range of values of 
that parameter that occur within the region defined by 
X 2 = Xmin + X 2 that we obtain when all parameters, ex- 
cept a, are allowed to vary. Here xt = 3-35 is the value 
of x 2 per degree of freedom for the best fit. This x 2 con- 
tour is that which yields la errors when the uncertainties 
in the individual frequency measurements are scaled to 
give x 2 = 1 ; assuming correlations in the errors of dif- 
ferent measurements are negligible. These estimates are 
given in column 4 of Table 2. These uncertainties do not 
reflect errors that could occur because (1) there could be 
correlations in the errors in the frequency measurements, 
(2) the flux model may not be precise, or (3) the model 
may not be formulated properly. We have attempted to 
crudely estimate the effects of the second possibility and 
the issue of possible systematic errors in general by vary- 
ing the parameters describing the flux model one at a 
time while fixing all other flux model parameters at the 
values given in Table 1 and then redoing the fit for each 
new set of values. In this exercise we did not limit the 
range of values of the modified parameter to those values 
that give a good fit to the flux data. The magnitudes of 
the variation of each derived frequency model parameter 
in this exercise are given in column 5 of Table 2. 

The fitted eccentricity of 0.033 indicates that the orbit 
is almost circular. Even so, a circular orbit fit (e = 0) 
results in a significantly poorer fit [xt = 4.57). The 
formal probability that the fit would yield e > 0.033 
when the true orbit is perfectly circular is ~ 1.1 x 10~ 7 ; 
this is equivalent to a detection of non-zero e with 5a 
confidence (cf. Lucy & Sweeney 1971). However, we 
note that the systematic uncertainty estimated in the 
exercise described above suggests that the significance of 
the measured eccentricity is closer to ~ 2.5cr. 

4.3. Pulse Arrival Time Model 

The pulse frequency fit takes advantage of the long 
baseline spanned by the observations, and yields an ac- 
curate determination of P rb- However, that analysis did 
not take advantage of the relative phases of the pulsations 
from observation to observation. We therefore undertook 
a pulse arrival time analysis of the PCA data obtained 
in the first epoch of RXTE observations. 

We first constructed a pulse profile for each observa- 
tion by folding data in 0.25-s time bins using the pulse 
period predicted by a preliminary but reasonably accu- 
rate frequency model that included both intrinsic fre- 
quency changes and orbital Doppler shifts. Profiles for 
selected observations are shown in Fig. 0] The pulse 
shape evolved as the flux increased and decreased. At 
lowest flux the profile exhibited a single peak, with a 
narrow dip spanning 0.2 in phase forming the primary 
minimum. As the flux increased, a shallower secondary 
dip around 0.25 later in phase developed, so that the 
profile exhibited one narrow and one broad peak. The 
primary minimum also became progressively broader and 
asymmetric, with a more gradual ingress and egress. As 
the flux decreased following the peak, the same system- 
atic variation in pulse profile was retraced in the opposite 
direction. 

For the determination of pulse arrival times, we folded 
data from 512-s segments. We defined the pulse arrival 
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time as the epoch of maximum flux, which was deter- 
mined by fitting a sinusoid to the pulse profile. The 
errors in the arrival times were estimated by varying the 
epoch of maximum flux in the fit over a grid to determine 
the interval in which A\ 2 < 1- 

To make effective use of these arrival times, one must 
determine the number of pulsation cycles between each 
measured pulse arrival time. Thus, we used the pre- 
liminary frequency model noted above to estimate cy- 
cle counts. When observations were separated by longer 
than ~ 2 d it was not possible to unambiguously de- 
termine the cycle count using the frequency data alone. 
Instead, we performed a joint fit of both the frequency 
and arrival time data. For each long data gap, we fitted 
both frequencies and arrival times before and after the 
gap for times as far forward and backward as possible 
without introducing additional cycle-count ambiguities. 
Using this approach, we were able to determine the cycle 
counts from MJD 51914 to MJD 52062. 

When the expression on the right in Eqn. 3 is substi- 
tuted into Eqn. [21and integrated, one obtains an express- 
sion for the pulse phase. Pulse arrival times may then 
be taken to be the times when the phase is an integer 
number of cycles. A constant of integration appears in 
going from frequency to phase; it may be taken as the 
arrival time of a given pulse. 

We then fit the arrival times within the range over 
which we can unambiguously phase-connect the obser- 
vations, as well as the frequency data before and af- 
ter. The best-fit model did not fit the data well, as 
indicated by the large value of reduced x 2 . One ma- 
jor reason for the lack of agreement between the model 
and the measured arrival times was an instance, around 
MJD 51956.4, where the pulse frequency changed rather 
rapidly (see Fig. |3J. The event is apparent in both 
the frequency measurements and in the arrival time (or, 
equivalently, pulse phase) residuals from a fit calculated 
using only pulse arrival times prior to the event. Ex- 
amination of pulse profiles revealed no significant pulse 
shape change around the time of the event. 

The residuals were consistent with a stepwise change 
of frequency of Af « +1.8 /iHz occurring around 
MJD 51956.41 (such a change is included in the frequency 
fit; see H4.2JI . In order to obtain limits on the magnitude 
and duration of this event, we fit a model to the arrival 
time data that included a step in frequency A/ g i that 
occurred via a constant frequency derivative over a time 
interval At g ] centered at time t g \. These parameters, as 
well as low-order polynomial terms in time, were allowed 
to vary in the fitting process while the orbital parameters 
and the torque coefficient were fixed at the values deter- 
mined from the best-fit frequency model. When this fit 
was performed on the arrival time data from the 7.8 day 
time interval MJD 51952.03 through 51959.83, we obtain 
A/ g i = 2.01 ±0.04 ^Hz, igj = MJD 51956.393±0.01, and 
Atgi < 0.4 days, where the errors have been scaled by the 
square root of \1 = 0.69. When this fit was performed 
on the arrival time data from a shorter time interval, i.e., 
MJD 51955.6 through MJD 51959.1, the best fit is ob- 
tained for A/ g i = 1.83±0.09 AiHz, t gl = MJD 51956.416± 
0.014, and At g i < 0.5 days, again where the errors have 
been scaled by the square root of xl = 0.67. The best- 
fit values of the parameters which describe the glitch 
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Fig. 5. — Pulse frequencies, phase residuals, and X-ray intensi- 
ties around the time of the glitch, (a) Measured pulse frequency 
corrected for the binary orbit with error bars corresponding to 3<r 
confidence intervals. (6) Phase residuals calculated from a best-fit 
model that does not include a component representing the glitch, 
(c) Phase residuals from a full model that does include a compo- 
nent representing the glitch. Note the difference in scale between 
panels b and c. (d) The 2-60 keV band source flux. 



do not appear to depend, within limits, on the length of 
the data set. The full sets of model parameters of both 
fits may be found in Table 3. We note that the parame- 
ters listed in Table 3 are consistent with an instantaneous 
change in frequency. 

5. DISCUSSION 

We measured the flux and pulse frequency evolution of 
the X-ray pulsar KS 1947+300 through a large outburst 
in 2000-01, as well as in a smaller outburst in 2002. We 
fitted the measured pulse frequencies and arrival times 
to derive orbital parameters, and to determine the tem- 
poral variation of the intrinsic spin frequency of the neu- 
tron star. We found that the pulsar revolves in a nearly 
circular orbit (e = 0.033 ± 0.013) with orbital period 
40.415+0.010 d, and projected semi-major axis 137+3 lt- 
sec. This orbital period is significantly shorter than 
that derived from ana lysis of the pre-outburst ASM light 
curve, 41.7 + 0.1 davs ijLevine fc Corbetll2000|) . Although 
it is clear that the value based on pulse timing is correct, 
the basis of the discrepancy is not clear. We note that the 
X-ray light curves of other Be-star/neutron star binaries 
sometimes do contain quasiperiodic variations that are 
asynchronous with the orbits (e.g., GS 0834-430, Wil- 
son et al. 1997; XTE J1946+274, Wilson et al. 2003). 

If we assume that the neutron star mass Mx = 1.4M©, 
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then the X-ray mass function of fx =1.71 ± 0.04 M & 
implies a companion mass of at least 3.4 M (for an in- 
clination i < 90°). We do not find eclipses in either the 
ASM or the PCA light curves, and can confidently con- 
clude that the system is not eclipsing, thus ruling out 
inclinations > 85°. Companion masses M c = 5, 10, or 
20 M , imply inclinations of 57°, 38°, or 28° respec- 
tively. This range of masses and inclinations is compati- 
ble with the BOVe spectr al type of the optical counterpart 
l|Negueruela et al.ll200^) . 

The wide orbital separation and low eccentricity place 
KS 1947+300 within the small class of neutron stars with 
high-mass binary companions which almost certainly did 
not receive large "kick s" during the sup ernovae which 
led to their formation ( Pfa hl et all [2002). The optical 
counterpart appears to be a Be star that is not highly 
evolved. Thus its radius is most likely much smaller than 
the radius of its Roche lobe, i.e., than 104, 141, or 191 lt-s 
for M x = 1.4 M Q and M c = 5, 10, or 20 M , respectively. 
Furthermore, the wide separation and relatively small 
size of the companion suggest that it is unlikely that tidal 
effects have circularized the orbit since the supernova 
event wherein the neutron star was formed. If the pre- 
supernova orbit was circular, the measured eccentricity 
e ~ 0.03 suggests that the mass ejected in the supernova 
must have been rather small, i.e., AM/Mfi na i < 0.03, 
where -Mfi na i, the post-supernova system mass, could be 
as large as 20M Q . The low eccentricity of this system 
clearly presents a number of difficulties. 

An outburst of KS 1947+300 was seen with BATSE 
in 1994, and no other outburst was reported until 
that of 2000-2001. BATSE operated continuously un- 
til 2000 May, while the RXTE ASM has operated more 
or less continuously since early 1996. BATSE would 
most likely have detected any pulsed flux greater than 
« 20 mCrab in the 30-75 keV band (equivalent to a flux 
2 x 10 -10 ergscm~ 2 s" 1 ) that was present for more than 
a few days. The ASM would most likely have detected 
any outburst longer than w5d that reached a strength 
of ss 1 ASM count s _1 (equivalent to 13 mCrab or to a 
2-12 keV band flux 3 x IP" 10 ergscm" 2 s" 1 )- Thus, as 
was also discussed bv iNeeueruela et all ll2003j) . we can 
rule out the occurrence of a moderate to strong outburst 
in the intervening time, although it is possible that weak 
outbursts could have been missed. 

From the end of the outburst in 1994 to the beginning 
of the outburst in late 2000, the pulsar frequency de- 
creased from 53.47 to 53.30 mHz. The average frequency 
derivative during this interval was / = —0.8 x 10~ 12 
Hz s _1 . Spin-down in quiescence likely arises from 
the dominance of magnetic torques contributed by field 
lines which thre ad the accretion disk out side the corota- 
tion radius fe.g. lGhosh fe Lamblll979alrjh . Alternatively, 
negative torques may arise from the "propeller effect" 
(111 arionov fc Sunvaevlll975(l . We note that the magni- 
tude of the quiescent spin-down was comparable to that 
measured for at least one ot her pulsar, A 0535 +26, of 
-(2.2 + 0.6) x 10~ 13 Hzs" 1 i|Finger et al.l[l99^l . While 
the values of bj for each epoch (see Table 2) are compa- 
rable to, or significantly larger than the quiescent spin- 
down, we caution that it is not correct to interpret these 
parameters independently from the torque model. The bj 
are coefficients of terms we used to augment the accretion 



torque integral in equation so as to more accurately 
model the intrinsic spin changes of the neutron star. The 
integral alone would give zero torque only at zero flux, 
whereas the various negative fitted values of the bj indi- 
cate that zero torque corresponds to non-zero flux at a 
different level for each epoch. 

The 2-60 keV flux of KS 1947+300 peaked at 9 x 
10~ 9 ergs cm -2 s -1 in 2001 February. This flux implies 
an X-ray luminosity of 1.1 x 10 38 ergss -1 for a source 
distance of 10 kpc. Given this accretion luminosity, a 
neutron star of mass Mx = 1.4 M®, radius R — 10 km, 
magnetic moment \i = 1 X 10 30 G cm 3 (equivalent to 
an equatorial surface magnetic field strength of 10 12 G), 
and moment of inertia J = 1 X 10 45 g cm 2 , a standard 
theoretical treatment (Eqn. 15 of lGhosh fc Lamblll979U: 
Doss fc Rappaportl 11976]) predicts / = 1.5 x 10" 11 Hz 
s _1 whereas the intrinsic frequency derivative around the 
time of maximum flux obtained from our pulse timing re- 
sults is / ~ 2.3 x 10~ n Hz s -1 . These values are consis- 
tent given the simple nature of the theoretical estimate, 
the possibility of deviations from the assumed neutron 
star parameters, and the uncertainty in the distance to 
the system. 

We found evidence for a particularly rapid variation of 
the pulse frequency around MJD 51956.41. The phase 
residuals around this episode were minimized by includ- 
ing a change in frequency of Af « +1.8 ^Hz that oc- 
curred at a uniform rate over a time interval At < 10 hr. 
The frequency change corresponds to a fractional change 
Af / f ~ 3.7 x 10~ 5 . If the frequency change occurred 
at a constant rate over a 10 hr interval, then the in- 
trinsic neutron star pulse frequency derivative must have 
reached at least 4.8 x 10 -11 Hzs -1 , more than twice as 
large as the maximum inferred from the torque model, 
2.3 x 10 -11 Hzs" 1 . We note that this event occurred 
close to the time of peak flux and, therefore, when the 
accretion torque was near its maximum value. While it is 
possible that the abrupt increase in frequency arose from 
a brief episode of enhanced accretion, we found no evi- 
dence for the increase in flux which might be expected as 
an immediate consequence (see Fig. |SJ). A transient flux 
increase could have been missed during the « 4 hr gap 
(MJD 51956.39 to 51956.56) between the observations 
closest to the estimated time of the glitch. However, the 
increase in X-ray luminosity required to explain the mag- 
nitude of the glitch would be a factor of S3 2, sufficient to 
exceed the Eddingon limit for a canonical neutron star 
at 10 kpc distance. Furthermore, four ASM dwells at 
MJD 51956.51 indicate a flux of 7.3 ± 0.6 count s _1 , con- 
sistent with the level before and after the gap. 

This event may have arisen from a sudden change in 
the spin frequency of the neutron star, analogous to 
the "glitches" observed in other types of neutron stars. 
However, there are also some significant differences be- 
tween the glitch in KS 1947+300 and those in radio 
pulsars and AXPs. First, in contrast to radio pulsar 
glitches, which frequently exhibit a change in the spin- 
down rate / following a glitch, we measured consistent 
values of / before and after MJD 51956.41. Second, 
the fractional change in frequency of 3.7 x 10~ 5 in the 
KS 1947+300 glitch is around an order of magnitude 
larger th an the typical value fo r "giant" glitches in radio 
pulsars ijKrawczvk et alJl2003|) . Finally, what is more 
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puzzling is that the frequency change in KS 1947+300 
was in the same sense as the mean spin-up around that 
time, i.e. the glitch resulted in a brief acceleration of the 
spin- up rate by at least a factor of 2 to 3. Since the neu- 
tron star was spinning up at the time of the glitch, one 
might expect that the core would be rotating more slowly 
than the crust, so that a glitch would cause a sudden de- 
crease in the pulse frequency. However, a long history of 
spin up and spin down episodes may have left the core 
spinning faster than the crust at the time of the glitch. 
The BATSE observations indeed show that the pulse fre- 
quency in 1994 was higher than the pulse frequency at 
the time of the glitch. 

This is the first reported glitch in an accretion-powered 
pulsar. If such events occur in other accreting neutron 
stars, why have they not been detected previously? The 
main reason may be that such events are simply difficult 
to detect; indeed, there are a number of conditions which 
are required for the detection of glitches in either accret- 
ing or rotation-powered neutron stars. The pulsars must 
be monitored over a period of days, with observations be- 
ing done at least several times per day. This is often more 
difficult to accomplish from space-based X-ray satellites 
than from ground-based radio observatories. Also, the 
extraction of pulse frequencies from previous X-ray ob- 
servations of pulsars have sometimes been done using 1 
day or longer integration time intervals, which makes it 
difficult to detect /iHz frequency changes occurring on 
time scales of hours or less. Sudden frequency changes 
could be lost among the other frequency changes from 
torque variations which typically take place on longer 
time scales. 

We can estimate the sensitivity for the detection of 
a sudden frequency change from a set of observations. 
For this purpose we assume that two estimates of the 
pulse frequency are made, one before and one after the 
time of a possible sudden change. For simplicity, we fur- 
ther assume that each frequency estimate is made from 
a set of N equally-spaced observations made over an in- 
terval of length T, that each observation yields a phase 
estimate with la uncertainty o^, and that there are no 
pulse counting ambiguities left unresolved. In this case 
the la uncertainty in the difference of the two frequency 
estimates is given by 
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where a is a factor of order unity. The length of the 
time interval T over which the frequency estimates are 
to be made is limited by the frequency changes generally 
ascribed to accretion torque variability. Equation^] indi- 
cates that the detection of a ~ 1 /iHz glitch requires, for 
example, N ~ 10 observations over the course of a few 
days with each observation yielding a modest accuracy 
phase measurement. 

The BATSE instrument on the Compton Gamma-Ray 
Observatory monitored many accretion-powered pulsars 
over the course of 9 years with a sensitivity of ~ 1 /xHz 
per day in favorable cases (Bildsten et al. 1997). Fre- 
quencies were reported for every day or every few days 
that pulsations were detected (Chakrabarty et al. 1997a, 
Finger et al. 1999, Koh et al. 1997, Chakrabarty et al. 
1997b). Though the BATSE results have not led to the 
report of any glitch, it is possible that some of the mon- 
itored pulsars did exhibit glitches that were not easily 
noticed because of the large frequency variations on time 
scales of days and weeks. 

A second possibility is that such events may also be 
quite rare. In this regard, we note that KS 1947+300 
was accreting near the Eddington limit when this glitch 
occurred. If glitches in accretion-powered pulsars only 
occur under such circumstances, then the opportunities 
for their detection would be quite limited, since the times 
when Galactic pulsars are observed to be accreting near 
the Eddington limit are relatively infrequent. Clearly, 
new analyses of archival data and additional observa- 
tions are called for in order to find other examples of 
glitches and to determine how often such events occur in 
accreting neutron stars. 
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TABLE 1 

Gaussian Components of the X-ray Light Curve 



Data set 


Component 
number 


Center (MJD) 


Width a (d) 


Peak Intensity 
(10 -9 ergscm _2 s _1 ) 


BATSE 


1 


49465.0 


10.0 


1.92 


RXTE 


2 


51855.0 


10.0 


1.92 




3 


51904.2 


10.0 


2.72 




4 


51945.4 


23.6 


6.43 




5 


51974.9 


32.1 


4.12 




6 


52020.5 


7.4 


0.449 




7 


52405.0 


17.0 


2.88 



a Standard deviation 



TABLE 2 

Frequency model parameters for KS 1947+30 



Parameter 


Epoch 


Fit Value 


Uncertainty a 


Uncertainty due 
to Flux Model b 


t i (MJD) 


BATSE 


49448.5 






/oi (mHz) 


BATSE 


53.4610 


0.0015 


0.0010 


6i (10- 12 Hz s- 1 ) 


BATSE 


-1.9 


0.9 


6.2 


t 02 (MJD) 


RXTE 1 


51869.397 






/ 02 (mHz) 


RXTE 1 


53.3000 


0.0003 


0.0005 


b 2 (10- 12 Hz s- 1 ) 


RXTE 1 


-3.4 


0.2 


0.6 


*03 (MJD) 


RXTE 2 


52362.115 






/03 (mHz) 


RXTE 2 


53.4448 


0.0008 


0.0036 


6 3 (10- 12 Hz s- 1 ) 


RXTE 2 


-7.2 


0.3 


3.5 


t g i (MJD) 


All 


51956.416 






A/ g i (Hz) 


All 


1.83 x 10~ 6 






a 


All 


0.75 


0.02 




P c 


All 


4.35 x 10~ 7 


0.07 x 10~ 7 


0.21 x 10~ 7 


PoA (d) 


All 


40.415 


0.007 


0.010 


Ox sin i (lt-s) 


All 


137.4 


1.2 


3.0 


T^/a (MJD) 


All 


51985.31 


0.07 


0.06 




All 


0.034 


0.007 


0.013 


I n 


All 


33 


3 




MM) (M ) 


All 


1.71 


0.04 





X 2 /A f (DoF) 329/98 



a The range of variation of the parameter values that occur within the \ 2 = Xmin ^min/^P°' 7 ) region that applies when all other 
frequency model parameters, except for a, are allowed to vary. 

b The range of variation of the parameter for the best fit when the parameters describing the flux model (Table 1) are varied one at a time 
(see text). 

c Torque coefficient in units of Hz day -1 



TABLE 3 

Pulse Arrival Time Models Around the Glitch 21 



Parameter Valuc b - Fit 1 Valuc b - Fit 2 

Time interval (MJD) 51952.03-51959.83 51955.6-51959.1 

to (MJD) 51869.397 51869.397 

fo (mHz) 53.29791 ± 0.00027 53.2976 ± 0.0014 

b (10~ 12 Hz s" 1 ) -15.6 ±0.2 -15.4 ±1.0 

tgi (MJD) 51956.393 ± 0.010 51956.416 ± 0.014 

At ! (days) < 0.4 < 0.5 

A/gi (^Hz) 2.01 ±0.04 1.83 ±0.09 

X 2 /N(DoF) 154.0/224 121.7/183 



a The following parameters were held fixed at the indicated values determined from the frequency model fitting: P or t, = 40.415 d; 
a x sini = 137.3 (lt-s); T^ /2 = {MJD) 51985.31; g = 0.018928; h = 0.027032; a = 0.75; fi = 4.35 X lO" 7 Hz day- 1 . The values of g and h 

correspond to an eccentricity e = y^(o 2 + h 2 ) = 0.033. The uninteresting parameter <f>o is not listed here. 

b Parameters without uncertainties were held fixed at the given values. The indicated uncertainties correspond to la confidence levels 
obtained from evaluation of A\ 2 = xi- Upper limits correspond to 2cr confidence limits. 



